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Abstract 

We study the impact of explicit chiral symmetry breaking of Wilson fermions on mesonic cor- 
relators in the e-regime using Wilson chiral perturbation theory (WChPT). We generalize the 
e-expansion of continuum ChPT to nonzero lattice spacings for various quark mass regimes. It 
turns out that the corrections due to a nonzero lattice spacing are highly suppressed for typical 
quark masses of the order aAQ CD . The lattice spacing effects become more pronounced for smaller 
quark masses and lead to non-trivial corrections of the continuum ChPT results at next-to-leading 
order. We compute these corrections for the standard current and density correlation functions. A 
fit to lattice data shows that these corrections are small, as expected. 
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I. INTRODUCTION 



The e-regime of QCD [H [2] offers various advantages for the numerical determination 
of the low-energy couplings (LECs) in chiral perturbation theory (ChPT), the low-energy 
effective theory of QCD. Only the leading order LECs, the pseudo scalar decay constant in 
the chiral limit F and the chiral condensate £, enter the predictions of ChPT through next- 
to-leading order (NLO) in the epsilon expansion. The Gasser-Leutwyler coefficients [2111] first 
appear at one order higher, thus making the e-regime attractive for precise determinations 
of F and E. 1 In addition, gauge field topology plays an important role in the e-regime [6]. 
ChPT makes predictions for correlation functions restricted to individual topological sectors, 
thus enlarging the number of observables that can be compared to numerical lattice QCD 
results. 

The role of topology and spontaneous chiral symmetry breaking has led to the widespread 
conviction that overlap [7J or domain-wall fermions [SJ [9J [10] are the preferred or even 
mandatory choice for lattice simulations in the e-regime. Consequently, a fairly large number 
of quenched simulations with these fermions in the e-regime can be found in the literature 
[TT| [121 [T3] HU [151 [H3 El UHl [19]. Even though the results were to a large extend promising, 
the main hurdle for progress in real QCD is the need for simulations with dynamical sea 
quarks, and these are extremely time-consuming. 2 So far only the JLQCD collaboration [21] 
has carried out a large scale dynamical simulation with overlap fermions in the e-regime, 
and the computer resources that went into this simulation are enormous. 

In contrast, fairly inexpensive simulations with tree-level improved Wilson fermions have 
been reported recently [22]. Reweighting as described in Ref. [23] has been used to reach 
small enough quark masses in order to be in the e-regime. The size of the box was L ~ 2.8fm, 
much larger than in all the simulations mentioned before. Quite surprisingly, the data 
for the axial vector and pseudo scalar correlation functions are very well described by the 
corresponding ChPT predictions, although chiral symmetry is explicitly broken for Wilson 
fermions. 

A similar observation has been made before by the ETM collaboration [2~i] |2"S] . Their 
data 3 , obtained with a twisted mass term [271 EH], also suggests that the e-regime can 
be reached with Wilson fermions. One might argue that automatic 0(a) improvement at 
maximal twist [29, 30, 31J may suppress chiral symmetry breaking effects. Still, the data 
obtained with Wilson fermions raises the question if and how the results can be interpreted 
in the presence of the explicit chiral symmetry breaking by the Wilson term. 

In this paper we study this question with Wilson ChPT [321 EH], the low-energy effective 
theory for lattice QCD with Wilson fermions. And indeed, our analysis suggests a natural 
answer to the question raised above. It turns out that the lattice spacing corrections are 
in general highly suppressed and show up at higher order in the epsilon expansion. For 



1 For a recent review of the various estimates see [5] . 

2 For a recent review see [20] . 

3 A review of the ETM results can be found in Ref. [26] , 



2 



example, for quark masses m ~ aAQ CD the deviations from the continuum results due to 
the 0(a) corrections enter first at next-to- next-to-leading order (NNLO). This is in contrast 
to the expansion in the p-regime, where the corrections already appear at next-to-leading 
order. This suppression is completely analogous to the suppression of terms involving the 
Gasser-Leutwyler coefficients in continuum ChPT in the e-regime. 

However, the power counting in WChPT can be different, depending on the relative size 
of the quark mass m and the lattice spacing a. For this reason we also consider a different 
power counting where m ~ aAq CD is no longer appropriate. In this case the lattice spacing 
corrections enter already at NLO. Interestingly, these corrections are entirely caused by the 
0(a 2 ) term in the chiral effective Lagrangian that also determines the phase diagram of the 
lattice theory [32] • The corrections linear in the lattice spacing, stemming from the effective 
Lagrangian and the effective operators, are still of higher order in the epsilon expansion. 
This is relevant in practice: The Wilson ChPT expressions contain only one more unknown 
LEC at this order, and the predictive power is not spoiled by a plethora of free fit parameters. 

We also use our WChPT results derived here for an analysis of the data of Ref. [22J. 
The corrections due to the nonzero lattice spacing turn out to be very small, supporting 
our theoretical analysis that these corrections are in general highly suppressed. This is 
very encouraging for lattice simulations with Wilson fermions. The impact of explicit chiral 
symmetry breaking for e-regime simulations is much less severe than previously thought. 
This makes simulations with Wilson fermions a serious and efficient alternative to those 
with chiral fermions. 



II. WILSON CHIRAL PERTURBATION THEORY (WCHPT) 
A. Chiral Lagrangian 

The chiral effective Lagrangian of WChPT is expanded in powers of (small) pion momenta 
p 2 , quark masses m and the lattice spacing a. Based on the symmetries of the underlying 
Symanzik action [SUES] the chiral Lagrangian including all terms of 0(p A ,p 2 m, m 2 ,p 2 a, ma) 
is given in Ref. [33]. The 0(ct 2 ) contributions are constructed in Ref. [36] and, independently, 
in Ref. [37] for the two-flavor case. 

In the following we will restrict ourselves to Nf = 2 with degenerate quark mass m. The 
continuum part of the chiral Lagrangian is the well-known Gasser-Leutwyler Lagrangian 
[21 [38] . With our notations the leading part reads (in Euclidean space-time) 

C 2 = ^Tr (dpUdprf) - ^mTr (U + tf) . (1) 
The field U containing the pion fields is defined as usual, 

U(x) = exp 0k(*)) , COO = C(x)T a . (2) 

The SU(2) generators are normalized such that Tr(T a T 6 ) = 5 ab /2, so T a = a a /2 in terms 
of the standard Pauli matrices a a . The coefficients B and F are the familiar leading order 
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(LO) low-energy coefficients. 4 Higher order terms are collected in the next-to-leading order 
Lagrangian £4 [3], which we do not need in this work. 
The terms involving the lattice spacing are as follows: 5 

C a = aW 45 Tr [dJJdprf) Tr (U + E7+) - amW^iTi (U + U ] )f , (3) 
= ^c 2 a 2 (Tr(U + U^)) 2 . (4) 

W / 45)W / 68 are LECs, similar to Gasser-Leutwyler coefficients in £ 4 . The quark mass and 
lattice spacing enter through the combinations 

fa = 2Bm, a = 2W a , (5) 

where Wo is another low-energy coefficient [33]. Its presence here is for dimensional reasons: 
Wo is of dimension three and a therefore of dimension two. Hence, a has the same dimension 
as the familiar combination Bm. 

We have chosen to parametrize the 0(a 2 ) contribution in terms of the LEC 0%. This 
coefficient plays a prominent role since its sign determines the phase diagram of the theory 
[52] . We briefly come back to this after we have discussed the power counting in WChPT 
in section III CI 

Note that the mass parameter m in eq. ([IJ is the so-called shifted mass [32] . Besides 
the dominant additive mass renormalization proportional to 1 /a it also contains the leading 
correction of 0(a). Consequently, the term F 2 aTr (U + W) /4 is not explicitly present in 
the chiral Lagrangian since it is absorbed in the shifted mass [33]. Of course, physical results 
expressed by observables do not depend on what mass is used for the parametrization of the 
chiral Lagrangian. 



B. Currents and densities 

The expressions for the currents and densities in continuum ChPT are well known. For 
example, the LO expressions for the axial vector current and the pseudo scalar density read 

= *Y Tr {T a (U^d,U - Ud^)) , (6) 

P c a t = i^Tr {T a (U - W)) • (7) 

In WChPT these expressions receive corrections proportional to the lattice spacing. The 
currents and densities in WChPT can be constructed by a standard spurion analysis, sim- 
ilarly to the construction of the chiral Lagrangian. One first writes down the most general 
current /density that is compatible with the symmetries using the chiral field U, its deriva- 
tives and the spurion fields. Then, in a second step, one imposes the appropriate Ward 

4 With our normalization the pion decay constant in the chiral limit is F w 93MeV. 

5 We essentially adopt the notation of Refs. [31] and [ST] . 
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identities valid in the theory. For the vector and axial vector current this has been done in 
[30] . Alternatively one can introduce source terms for the currents and densities and con- 
structs the generating functional, as has been done in Ref. [32]. Carrying over the notation 
of this reference the axial vector current including the leading corrections of 0(a) reads 

^,wchPT = ^,cont (l + ^aW^Tt (U + U^+ 2aW w d,Tr (T a (U - 17+)) . (8) 

This axial vector does not satisfy any specific renormalization condition. Imposing a par- 
ticular renormalization condition leads to a finite renormalization. Explicitly, one introduces 

m 

^,re» = Z A A^ WCbpT ( X ) (9) 

with a renormalization factor Za- 6 Now one can impose a renormalization condition and 
demands that it is satisfied by A" n ; this determines Za- For instance, in Ref. [40J the 
massless chiral Ward identity has been imposed. 
Quite generally, Za has the form (up to 0(a)) 

16 

Z A = 1 + — 2 aW A (10) 

with an unknown coefficient Wa- This form reflects the fact that, by construction, the 
WChPT current reduces to the correct continuum current for a — > 0. Consequently, Za is 



equal to one in the continuum limit, and (10) is the leading generalization for a ^ 0. Hence, 



the general form of the renormalized current is, up to O(a), 

^,WChPT = ^,cont (l + ^ ^ 45Tr ( U + W ) + AWa \) 

+2aW w d,Tr(T a (U-U^)). (11) 

For brevity we have dropped the subscript "ren" on the left hand side. Terms of 0(am, a 2 ) 
will be present at higher order in the chiral expansion. 
Analogously, we use the expression [39] 

^WChPT = ^ont (l + + U ) + AWP \ ) (12) 

for the pseudo scalar density. The contribution involving the LEC Wp (not present in Ref. 
|39j) stems from a general renormalization factor Zp = 1 + 16aWp/F 2 , which we also allow, 
even though the results derived in this paper will not depend on the details of the O(a) 
correction. 

6 Za is a renormalization factor in the effective theory and should not be confused with Za in the underlying 
lattice theory. 
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C. Power counting in infinite volume 



In WChPT there are two parameters that break chiral symmetry explicitly, the quark 
mass m and the lattice spacing a. The power counting is determined by the relative size of 
these two parameters. 

The literature [391 l4"T] distinguishes two quark mass regimes with different power count- 
ings: (i) the GSM regime 7 where m ~ aAq CD and (ii) the Aoki regime where m ~ a 2 AQ CD . 
A priori one does not know in which regime one actually has performed a simulation. For 
this to decide one has to compare with the predictions of WChPT and check which expres- 
sions fit the data better. However, recalling how lattice simulations are typically done one 
can easily imagine that one starts in the GSM regime and by lowering the quark mass at 
fixed lattice spacing one will eventually enter the Aoki regime. 

Depending on the particular regime, the LO Lagrangian is different. Since m ~ a 2 Aq CD 
in the Aoki regime, also the C a 2 part in Q counts as LO [37] : 



The effects due to a nonzero lattice spacing are much more pronounced in the Aoki regime. 
Non-trivial phase transitions become relevant [32] and additional chiral logarithms propor- 
tional to a 2 appear at one loop [33, 112] • 

D. The pion mass and the PCAC mass in infinite volume 

It is useful to derive the pion mass and PCAC mass at LO. 
We start with the calculation of the pion mass. Expanding the LO chiral Lagrangian to 
quadratic order in the pion fields we obtain 



The sign of c 2 determines the phase diagram of the theory [32J. 8 For c 2 > there exists 
a second-order phase transition separating the Aoki phase [13]. The charged pions are 
massless in this phase due to the spontaneous breaking of the flavor symmetry. The pion 
mass vanishes at m = C2d 2 /B. For even smaller values of m the charged pions remain 
massless, while the neutral pion becomes massive again [32J. 

Negative values of c 2 , on the other hand, imply a first order phase transition with a 
minimal non-vanishing pion mass. All three pions are massive for all quark masses, and the 
pion mass assumes its minimal value at m — 0, resulting in 



GSM regime : 
Aoki regime : 




GSM regime : 
Aoki regime : 



M 2 = 2Bm, 

Ml = 2Bm-2c 2 a 2 . 



(13) 
(14) 



Ml 



0,min 



2\c 2 \a 2 . 



(15) 



7 GSM stands for generically small masses. 

8 Note that our definition for C2 differs by a factor F 2 a 2 from the one in Rcf. [32 . 
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Note that magnitude and the sign of C2 are a priori unknown and depend on the details 
of the underlying lattice theory, i.e. what particular lattice action has been used. 
The PCAC quark mass is defined by the ratio (no sum over a) 

(d,A^x)P a (0)) 
mpCAC " 2(P-(x)P-(0)> ' (16) 

where angled brackets indicate expectation values. Expanding the current and the density 



in Eqs. (11) and (12) to 0(£ a ) we find 



Afa) = -iFd^ a {x) (1 + ac A ) , (17) 

p a ( x ) = iFBC(x) (1 + ac P ) . (18) 

where here and in the following we for simplicity no longer write subscript "WChPT" . We 
also introduced the short hand notation 

16 

c A = ^2W [W 45 + W A - W 10 /A] , (19) 
16 

c P = j^2W [W 68 + Wp], (20) 
for the combinations of LECs in the effective current and density at this order. The corre- 



lation functions in (16) are now easily computed at LO, yielding 

M 2 ( \ 
mpcAC = 2B V + a( ° A ~ ° P V • (21) 

Using the tree-level pion mass obtained above we find 

GSM regime : ^pcac = m{\ + a{c A — Cp) \ , (22) 
Aoki regime : m PC AC — (m — ^ fl2 ) ^1 + a ( c A ~ Cp)j . (23) 

In the GSM regime the result is rather simple and the PCAC mass is equal to the shifted 
quark mass, up to corrections of 0(ma). This is no longer true in the Aoki regime. Still, 
the contribution proportional to (c A — cp) is subleading and can be ignored if one works to 
leading order in the quark mass. 



Equations (22) and (23) allow to replace the shifted mass m, which is just a parameter 
in the chiral Lagrangian, with the PCAC quark mass. The latter is an observable which is 
often used in lattice simulations. 

Note that the results above reproduce a well-known fact, namely that the PCAC mass 
depends on the particular renormalization conditions imposed on the axial vector current 
and the pseudo scalar density. Different renormalization conditions show up as different 
values of c A and cp. 
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III. WCHPT IN THE EPSILON REGIME 



A. Continuum ChPT in finite volume 

Consider continuum QCD with Nf degenerate quark masses in a hypercubic volume 
V = TL 3 , with T,L^> 1/Aq CD . Finite-size effects can be systematically studied by means 
of the corresponding chiral effective theory [TJ |2j HI] . In this section we summarize the main 
aspects of finite-volume chiral perturbation theory in the continuum. 

If the pion Compton wavelength is much smaller than the size of the box, M n L ^> 1, 
finite-volume effects can be treated in the chiral effective theory by adopting the standard 
p-expansion, where the power-counting in terms of the momentum p is given by 

m ~ 0(p 2 ), 1/L, 1/T , d, ~ 0(p), e ~ 0{p). (24) 

For asymptotically large volumes, one expects the finite-volume effects to be exponentially 
suppressed by factors ~ e~ M ^ L . 

On the other hand, approaching the chiral limit by keeping [i = mY,V < 0(1) (but still 
L ^> 1/Aqcd), where £ = F 2 B is the quark condensate in the chiral limit, one explores the 
domain where the pion wavelength is larger than the size of the box, M n L < 1. In this case 
the pion zero-mode gives a contribution to the propagator proportional to 1/MqV, which 
cannot be treated perturbatively but has to be computed exactly [HE]. This is achieved by 
factorizing the pseudo Nambu-Goldstone boson fields as 

U(x) = exp0j£Or)W (25) 

where the constant Uq 6 SU (Nf) represents the collective zero-mode. The nonzero modes 
parametrized by £, on the other hand, can still be treated perturbatively and satisfy the 
condition 

/ d 4 x£(x) = 0, (26) 
Jv 

since the constant mode has been separated. 

The zero-mode contribution proportional to 1/MqV in the pion propagator diverges in 
the chiral limit and a reordering of the perturbation series that sums all graphs with an 
arbitrary number of zero-mode propagators is necessary [2]. This reordering is achieved 
with the power counting 

m~0(e 4 ), 1/L, 1/T, d, ~ 0(e), £ ~ O(e). (27) 

Mass effects are suppressed compared to the p-regime, while volume effects are enhanced 
and become polynomial in (FL)~ 2 . Since Mq is proportional to m the combination 1/MqV 
now counts as e°. Consequently, all graphs that exclusively involve zero-mode propagators 
count as 0(1) and are unsuppressed. The key point here is, that the counting of the quark 
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mass is dictated by the counting of L by demanding 1/MqV = 0(e°). We will use this in 
the next section in order to establish the counting rules in WChPT. 

With the factorization given in Eq. (25), the leading order continuum partition function 
in the e-regime is given by 



SU(N f ) 



[,11 ■ ,] / [<| oxp { \ I d*xTr(d^d,0 + ^Tr(f/ + Ui)y (2S) 



The integration over the perturbative degrees of freedom [<i£] gives rise to the usual Wick 
contractions, while the zero-mode integrals over [dUo] must be computed exactly. Notice 
that by going to 0(e 2 ), by factoring out the constant zero-mode from the measure, one 
obtains 

[dU] = [d£][dU ] (1 + A(0 + 0(e 4 )) , (29) 

with 



_ 2N f 1 

for a general value of Nf [2j H5] . 



M0 = -^ V l v d 4 xTr(e(x)) (30) 



B. Power countings for the epsilon regime in WChPT 

Like the continuum effective theory, WChPT can also be formulated in a finite volume, 
in particular the e-regime discussed in this section. 

In WChPT we have additional low-energy constants and the lattice spacing as an additional 
expansion parameter. The main task is to decide how to count these in the epsilon expansion. 

Just as the continuum LECs F and S, we count all the additional LECs associated with 
the lattice spacing to be of order e°, 

c 2 ,c A ,c P ~ 0(1). (31) 

The counting of the lattice spacing a is more complicated. The general strategy is to follow 
the infinite-volume procedure and determine the power counting depending on the relative 
size of m and a. At finite volume, once the counting of m is fixed by the counting of L, we 
obtain the counting of a. 

We start with the GSM regime. The LO Lagrangian and the pion mass Mq are as 
in the continuum e-regime, so we conclude m ~ 0(e 4 ) by the same arguments as in the 
previous section. Since the GSM regime is defined by m ~ aAQ CD we are immediately lead 
to a ~ 0(e 4 ). 

The Aoki regime is more subtle. According to our assumption, the pion mass Mq, given 



in (14), is now a sum of two terms of equal order. If c 2 < 0, it is a sum of two positive 
terms. Hence, a small pion mass of order 0(e 4 ) implies that both terms, 2Bm and 2\c 2 \c? 
are small too and also of order 0(e 4 ). 
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If C2 is positive, the pion mass is the difference of two positive contributions. This leaves 
the possibility that M5 2 is small, even though the individual terms 2Bm and 2\c2\a? may 
not be small and only their difference is. A pion mass of order e 4 may be obtained by the 
difference of two order e 2 or e 3 terms, for example. 

We do not think that this is a likely scenario. Present day lattice simulations are usually 
done with small lattice spacings less than 0.1 fm and the 0(a 2 ) corrections are expected to be 
small in this case. Hence, we assume that a 2 ~ 0(e 4 ) in the Aoki regime. This assumption, 
together with the requirement Mq ~ 0(e 4 ) then also leads tom~ 0(e 4 ), the same counting 
as for c 2 < 0. 

The epsilon expansion allows us to introduce yet another regime where we count a ~ e 3 . 
Just by the powers of e this is an intermediate regime between the GSM and Aoki regime. 
One may think about it as the GSM regime but at a larger lattice spacing (or smaller quark 
mass). Its usefulness will become clear in the next section when we discuss the epsilon 
expansion of correlation functions. 

All three countings we introduced are well defined and are appropriate for a particular 
relative size between m and a. In order to be able to refer to these regimes we introduce the 
following nomenclature: 

GSM regime : a ~ 0(e 4 ) , 

GSM* regime : a ~ 0(e 3 ) , (32) 
Aoki regime : a ~ 0(e 2 ) . 

For fixed values of m and a in a given regime, one can match for instance the time-dependence 
of current correlators with lattice QCD results in order to extract the corresponding LECs. 



C. Epsilon expansion of correlation functions 

We will be interested in correlators of the pseudo scalar density and the axial vector 
current. These correlators have been calculated before through NNLO in continuum ChPT 
[35] . In powers of e this corresponds to 0(e 4 ) for the {P a (x)P a (0)) correlator and 0(e 8 ) for 
(Afa)A*(0)). 

In order to discuss the epsilon expansion in WChPT let us split an arbitrary operator 
and the action in WChPT into the continuum part and a remainder proportional to powers 
of a, 

O(x) = O ct (x) + SO(x) , (33) 
S = S ct + 5S. (34) 

Expectation values are generically defined as 

(O) = i f [dU}e~ s O, (35) 
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where Z is the partition function 



[dU]e- s . (36) 

The two-point correlator (0\{x)02(fy) = (O1O2) (for notational simplicity we suppress the 
dependence on x) can then be written according to 

(Ox0 2 ) w = (Oi,ct02,ct> + 5(Oi0 2 >, (37) 
5(0x02} = (O hct 50 2 + 80t0 2 , c% ) - (0 ltCt 2 , ct 5S) + (0 1>ct 2>ct )(5S) . (38) 

Here we have approximated exp(— SS) ~ 1 — 5S and we dropped all higher corrections. Note 



that the expectation value on the left hand side of (37), labelled with a subscript "W", is 
defined with the full action S in the Boltzmann factor, while on the right hand side it is 
defined with S ct only (for notational simplicity we suppress a subscript "ct"). 

The discretization corrections for the pseudo scalar and the axial vector can be read 



off from (11) and (12). For what matters here we can simplify these expressions. We are 
interested in the power counting for the epsilon expansion, and for this we can ignore all 
constants which count as 0(1). Therefore, we write 

5P a oc a(Tr (U + U ] ) + l)P c a t . (39) 

SP a is proportional to the continuum density itself. As mentioned before, the epsilon expan- 
sion of P c " starts with 0(e°). The "scalar density" Tr(U + W) also starts at 0(e°). Hence, 
by considering the continuum contribution at LO, the leading correction in the epsilon ex- 
pansion due to lattice terms is completely determined by how we count the lattice spacing 
a. In the last section we defined three different countings, so for now we leave it unspecified, 
write a ~ e Ua where n a counts the epsilon powers for a, and obtain 

5P a ~ e n \ (40) 

Note that the symbol ~ stands here just for the leading lattice contribution in the epsilon 
expansion. 

Analogously, we find for the axial vector (dropping again irrelevant constants) 

5A« oc a[(Tr (U + f/t) + 1)^ + «9,Tr (T a (U - f/t)) ] . (41) 

Both, A£ )Ct and <9 M Tr {T a (U - f/t) have an open Lorentz index and, therefore, contain at 
least one derivative acting on at least one power of ^(x). Hence, their continuum epsilon 
expansion starts at 0(e 2 ) and we find for the leading lattice corrections 

8AI ~ t na+2 . (42) 

Finally, we have to take into account the lattice corrections due to the contribution 5S. It 
will be useful to split them into two parts. The first one, denoted by 5S a , contains the terms 
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of C a in ([3]). These terms start at 0(e n ") (having taken into account e 4 from the integration 
over space-time), 

SS a ~ e n ". (43) 

The second contribution, 5S a 2, contains only the a 2 term proportional to C2. Therefore, it 
counts as 

SS a2 ~ e 2 ^" 4 . (44) 

After these preparations we can determine at which order the lattice spacing effects enter 
the PP and the A A correlator. 



1. GSM regime 

In the GSM regime we set n a = 4 and find 

5(P a (x)P a (0)) ~ 0(e 4 ), (45) 
5(^(x)^(0)) ~ 0(e 8 ) . (46) 

The corrections due to the lattice spacing first affect both correlators at NNLO. Up to 
NLO the results obtained in continuum ChPT are the appropriate ones. This is quite 
remarkable and may explain why numerical data generated recently [22] could be fitted very 
well using the NLO continuum expressions. Note that this suppression to NNLO holds for 
the unimproved theory. The reason is that the terms linear in a are accompanied by at least 
one additional power of m or <9 M £, and are therefore of higher order. 



2. GSM* regime 

Here we have n a = 3 and obtain 

5(P a (x)P a (0)) ~ 0(e 2 ), (47) 
5(Al(x)A«(0)) ~ 0(e 6 ) , (48) 

hence the corrections enter at NLO. Interestingly, the dominant term here comes only from 
the correction 5S a 2. The other corrections from the 0(a) contributions in the currents, 
densities and SS a start at e 3 and e 7 , respectively. Therefore, they are of higher order in the 
epsilon expansion, even though they are still lower than the NNLO contributions, which 
start at e 4 and e 8 , respectively. 

Notice that in the GSM and GSM* regimes the leading order partition function is like the 



continuum one given in Eq. (28). In particular, the exact zero- mode integrals are computed 



with respect to the same Boltzmann factor as in continuum ChPT. 
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3. Aoki regime 



The modifications in the Aoki regime are more pronounced than in the previously dis- 
cussed regimes. Here, cut-off effects show up already at LO. Even worse, the corrections 
can no longer be linearly added to the continuum result. The reason is the correction SS a 2, 
which gives a zero-mode contribution of order e°. Hence, it is no longer justified to expand 
completely the exponential exp(— S a 2) « 1 — S a 2. The zero-mode contribution of order e° 
has to be included in the leading order Boltzmann factor. That is, the partition function 
becomes 



J Aoki 



f [dU ] [ [dfl exp ( \ f d 4 xTr(d^d,0 + ^Tr(C/ Q + J) (49) 

JSU(2) J V 1 JV 1 



c 2 F 2 a 2 V _ t rri ..r 



16 

This modification affects all constant integrals and probably leads to non-trivial changes of 
the continuum results. Note that the other 0(a) corrections (from 50 and 5S a ) are of order 
e 2 and show up at NLO only. 



D. Comment on O(a) improvement 

The results in the previous section are valid for unimproved Wilson fermions. It is natural 
to ask how (non-perturbative) 0(a)-improvement changes these results. 

If the theory is non-perturbatively improved the corrections 50 and 5S a are absent, 
and modifications are caused by 5S a 2 only. We have seen that this term is the dominant 
correction and the others are subleading. Consequently, the epsilon expansion is essentially 
unaltered for the improved theory, since only subleading terms vanish. 



IV. LEADING CORRECTION IN THE GSM* REGIME 

We already mentioned in the introduction that the epsilon expansion is advantageous for 
the determination of F and S, since the Gasser-Leutwyler coefficients first enter the ChPT 
formulae at NNLO. The same is true for WChPT in the GSM regime, where the additional 
lattice spacing contributions enter at NNLO too. In other words, working through NLO the 
results for the GSM regime are the same as those in continuum ChPT. 

The first non-trivial modification of the continuum NLO results appears in the GSM* 
regime, and in this section we compute the leading correction to the PP and AA correlator; 
the results for some other correlators are given in appendix [Cj This correction is caused by 
the constant term of the 5S a 2 contribution, 

6(0^)0^)) = -(0^(x)0^%(y)5S a 2) + (0^(x)0^ t (y))(5S a 2), (50) 

leading 
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where 



5S a 2 = f-(Tr(U + U )) 2 . (51) 
lo 

The superscript "LO" refers to leading order in the e-expansion and we have introduced the 
dimensionless quantity 

p = F 2 c 2 a 2 V. (52) 
Notice that if (0\° t (x)0%° t (y)5S a2 ) = {0^ t (x)0^ t (y))(5S a 2), i.e for disconnected inser- 



tions, the leading correction in Eq. (50) vanishes. This happens for instance for left handed 
(V-A) current correlators jlB] , and more general for correlators which do not get zero-mode 
contributions at leading order. 

The epsilon expansion in continuum ChPT allows predictions for fixed topological sectors 
[0] . Since chiral symmetry is explicitly broken in lattice QCD with Wilson fermions, an exact 
definition for the topological charge does not exist at nonzero lattice spacing. For this reason 
we will only give results for observables where the sum over all topological sectors has been 
performed. 



A. Preliminaries 



Calculations of mesonic 2-point functions in the e-regime have been pioneered in Ref. 
[4"5] . The integration over the non-constant modes £(x) is done perturbatively as in ordinary 
chiral perturbation theory in the p-regime. We summarize the corresponding propagators 
and other useful properties in Appendix [A} 

The integral over the constant mode Uq has to be done exactly. In our particular case 
with Nf = 2 we encounter integrals of the type 

(9(Uo)) = ±-f [dU ] g(U ) e^ u ° +u ^ , (53) 

^0 JSU{2) 

where Z$ is the continuum partition function associated to the zero modes, 



[dU }e^ Uo+u ^ , (54) 

SU{2) 

and fj, denotes the standard combination 

fi = mZV. (55) 



Quite generally, the integral (53) leads, at least for the types of g we are considering, to 
expressions involving modified Bessel functions I n (z) with integer index n. They satisfy 
numerous recursion relations [47j which allow us to express all integrals in terms of two 
Bessel functions, which we choose to be I 2 and I\. In Appendix [B] we collect various integrals 
that one encounters in calculating the PP and AA correlator in the GSM* regime. 
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As an example let us consider the expectation value of the quantity 5S a 2 } defined in Eq. 



(51), that is part of the correction in (50). By using the integrals given in the Appendix we 



obtain 



(SS a2 ) = p 1 



3 I, 



2 / z/ 1 (2 A o; • 



(56) 



B. The PP correlator 



For the PP correlator we introduce the definition 

(P a (x)P b (y)) = 5 ab C PP (x - y) 



(57) 



which takes into account translation invariance and the trivial dependence on the flavor 
indices. In the GSM* regime, Cpp(x — y) can be written through NLO as the sum of the 
continuum correlator and a correction proportional to a 2 , 



C PP (x - y) = C PP)Ct (x - y) + C PP ^{x - y) . 



(58) 



The continuum correlator for generic Nf at NLO (which corresponds to 0(e 2 )) is given by 
(see also [15] ) 

Cpp, ct (x - y) = C P + a P G(x - y), (59) 



where G(x — y) is the finite-volume massless scalar propagator defined in Eq. (Al) and 



C f 



ap 



Y 2 

^eff 



8(N] - 1) 
E 2 

4F 2 (7V 2 - 1} 



(Tr[(C/ - ^) 2 ])e ff - ^([Tr(f/ - W> 



N 



cff 



2JV?-4-— (Tr(C/ 2 ) + Tr(C/t 2 )) 



(60) 



2 , /iV 2 + l 



((TrU o y + (Tr^) 2 ) 



Nf 



Nf 



(61) 



The expectation values with the subscript "eff" are defined like in Eq. (53) with fi replaced 
by 

/i eff = mS eff V, (62) 
where S e fj is the quark condensate at one loop [45] 



^eff 



S 1 + 



Nf-l l Pi 

N f F 2 ^v)- 



(63) 



Pi is a so-called shape factor and is defined in Eq. (A5). 



For the particular case Nf = 2, after the explicit computation of the zero-mode integrals 
according to Appendix |Bj one gets 

X 2 



c, 



ap 



cff J 



2/ieff A(2/Ueff) 



2F 2 



2 - 



l / 2 (2m) " 
M/i(2m). 



(64) 
(65) 
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For the leading lattice correction to the continuum result, as given in Eq. (50), we find the 
0(e 2 ) contribution 

C PPt a2 = p—App, (66) 



with 



\pp 



5/i/ 2 (2/i) - 10/i(2/x)/ 2 (2/x) - 3/x/|(2/x) 



(67) 



Interestingly, the correction App is finite in the limit \i — > 0, as is easily checked using the 
leading order Taylor expansions for the Bessel functions jH], ii(2/x) ~ \x and l2{2p) ~ A i2 /2. 
We are not aware of a rigorous argument that App has to be regular at vanishing /1, since 
this correction ceases to be valid for small enough quark mass where one enters the Aoki 
regime. A singularity at // = would have been a clear signal for this breakdown of our 
calculation, however, this signal is not present in the result, at least not at the order in the 
chiral expansion we are working here. 

For the matching with numerical results obtained in lattice simulations we are interested 
in the correlation function integrated over the spatial components, 



C P p(t) = / d 3 xC PP (x - y) 



Cpp,ct(t) + 



3y2 



-pA 



pp 



where 



Cpp )C t{t) 



L 3 E 2 ff J 2 (2/i 



off J 



TS 5 



2 fi eS h(2ti eS ) 2F 2 



htit/T) 



2 - 



l/ 2 (2/i) 
///i(2/i)_ 



(68) 



(69) 



The time dependence is given by the parabolic function hi defined in Eq. (A9). 



The AA correlator 



The AA correlator is computed along the same lines. For simplicity we consider the 
time-component correlator and we define 



(A a (x)A b (y)) = 5 ab C AA (x - y) . 



(70) 



Similarly to the PP correlator we split Caa at NLO into a continuum part and a correction 
proportional to the lattice spacing, 



C A a{x - y) = C AAtCt (x -y) + C AA ^(x - y) . 



(71) 



The continuum contribution at 0(e 6 ) for x 7^ y and generic Nf has been calculated before 
(see also [49J) and is given by 



C AAtCt (x -y) = a A d xo d yo G(x - y) + P A K 00 (x - y) + ^ A d XQ d yo H(x - y), 



(72) 
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where the following short hand notation has been introduced: 



a A 

0A 
1A 



F 2 N f f3 1 



N f 



(2 - (Jo)) 



(nUo + Ut))^-. 



(73) 

(74) 
(75) 



The functions K^ v and H are given in Eqs. (A7) and (A8). Moreover, we have introduced 
the quantity 



Jo 



1 



TrC/oTrC/J + Nf - 2 



(76) 



Like for the PP correlator, the subscript "eff" refers to the substitution p — > p eS in the 
zero-mode integrals. For the particular case we are considering, Nf = 2, the results (73) - 
(75) reduce to 



a A 
Pa 

1A 



cS 



/^eff A(2/X e ff)_ 



+ 2 



01 



/x/i(2/x)J 



2 / 2 (2/x) 
2/i/ 2 (2/i) 



(77) 
(78) 
(79) 



J 1 (2//) ' 

In analogy to the PP correlator, the 0(a 2 ) contribution can be computed according to Eq. 
(50), and we obtain 



C A a a A x -y) = -^d xo d yo G{x - y)pA AA , 



with 



^AA 



-5/i/ 1 2 (2/i) + 10J 1 (2/i)J 2 (2/i) + 3/i/ 2 (2/i)- 



-2A 



pp 



(80) 



(81) 



Note that this correction affects only the coefficient a A in Eq. (77), which will be modified 
by the presence of lattice artifacts. 

The 0(a 2 ) correction is, up to a sign and a factor two, the same as the correction for 
the PP correlator. As far as we can see there is no deeper reason for this. It is simply a 
consequence of the fact that the zero-mode integrals for the leading order continuum PP 
and AA correlator are very similar, both lead to the same contribution involving the ratio 
J 2 (2//)//i/ 1 (2/i). 

By integrating over the spatial coordinates and using the properties listed in Appendix 
|A~1 we obtain for t ^ 



C AA {t) = / d 3 xC AA (x - y)\ y =o = C AA>ct (t) 



IT 



pA 



AA, 



(82) 
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where the continuum for N 



CAA,ct{t) 



+ 



2 explicitly reads 
T 



1 T 

--a A + — k 00 (3 A 



2T 



V 



1 



cS 



T 



1 



/ 2 (2/x) 
V *Vi(2/x) 



TVV 
2T ninA nh(2n) 



ta x {2n)) 



Here fcoo is another shape factor defined in the appendix, Eq. (A6). 



(83) 



D. The PCAC mass 

The correlators in the previous section are given as functions of m, the shifted mass. 
This is the mass parameter in the chiral Lagrangian and a priori not an observable. Here 



we compute the PCAC mass, defined in (16), and use it in the next section to replace m 
with mpcAc- 



We have already calculated the denominator of (16), and the numerator can be done 
analogously. Let us define 



{d,A^x)P\y)) = 5 ab C dAP (x-y) 

To leading order in the epsilon expansion we find the result 

S/ 2 (2/i) 



c d AP,ct( x - y) 



Vh(2fx) ' 



Dividing this by the leading order result of 2Cpp ct in Eq. (59) we obtain 



"^PCAC 



m . 



(84) 



(85) 



(86) 



This is just the result of continuum ChPT, where it is not surprising because the PCAC mass 
stems from the PCAC Ward identity. Note, however, that both numerator and denominator 
contain non-trivial Bessel functions which cancel in the ratio. This cancellation will no 
longer happen with the lattice spacing corrections included, since the PCAC relation no 
longer holds. 



The leading correction to the numerator in the GSM* regime is given by (50) with 0\ 



d fl A a (x) and O2 = P a (y) (no sum over a). The computation is straightforward as the ones 
in the previous sections and we find 



C dA p(x - y) 



S / 2 (2/x) 
Vh(2fx) 



3p_ 
2fi 2 



/// 1 (2/i) /x/ 2 (2//) \ 
I 2 (2n) 7 1 (2/i) J 



(87) 



Dividing by 2Cpp given in (58) we obtain the leading 0(a 2 ) corrections to the PCAC mass: 

1 + P 



T^PCAC 



m 



(88) 



Ji(2aQ V 

The key observation here is that the PCAC mass is equal to m, up to a correction of 
0(ma 2 V), which is e 2 higher in the epsilon expansion in the GSM* regime. 
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E. The correlators as a function of the PCAC mass 



The final step we have to do is to replace m by mpcAc i n the correlators. We first invert 

2 h(2p) V 



result (88) 



H = H 



1-p 



p 2 p,I 2 (2p,)J 



where 



H = rapcAcSV. 



(89) 



(90) 



In the NLO contributions of the correlators we can simply replace m = mpcAC; H — since 
the corrections are higher than this order. In the LO term, however, we have to use the full 



expression (89), which gives rise to additional corrections proportional to p. 



Eq. (89) has to be inserted into the Bessel functions I n (2p). Since the correction propor- 
tional to p is e 2 higher in the epsilon expansion we can Taylor-expand, 



J n (2/x) = / ft (2A)-2p/2/;(2/i) 



P 2 



him \ 
Phm) 



(91) 



and drop the higher order terms. The final results for the correlators can be brought into 
the form 



C PP {t) 
Caa if) 



L 3 £ 2 

C P p iCt (t) H — pA a 2 

F 2 

C A A,ct{t) + ^rpA a 2 , 



where the continuum correlators are as in (69) and (83), but with the replacements \i 
and /i e ff — ► fieS, with 

fieS — ^PCAcSeff^ • 



(92) 
(93) 



(94) 



The correction A a 2, which depends on jl, captures the lattice spacing artifacts and reads 
4/i 2 / 1 3 (2/i) - lljll 2 (2jl)l 2 (2jl) + 2(3 - 2p 2 )h(2p)I 2 (2p) + 5/LTf(2/i) 



2/23i?(2/i)/ 2 (2/2) 



(95) 



Note that it is regular at jl = 0. 

Eqs. (92) and (93) are our final results for the GSM* regime. (In Appendix [C| we also 
give the corresponding expression for the vector current correlator.) These results are re- 
markable and perhaps surprising in two ways: (i) The ft dependence of the 0(a 2 ) correction 
is identical for both correlation functions, (ii) Besides the continuum LECs F and S only 
one more unknown LEC enters these expression, the parameter c 2 - The second feature is 
very advantageous in practice when our results are used to fit numerical lattice data. 

A different question is the actual size of the 0(a 2 ) correction, which is directly propor- 
tional to c 2 . In the next section we try to give at least a rough answer to this question. 
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V. NUMERICAL TESTS 



A. General considerations 

For the pseudo scalar and axial vector correlators, the leading 0(a 2 ) correction in the 
GSM* regime is just a shift of the constant part. The question is how big this correction is 
in a typical e-regime simulation. As a measure for the correction we study the ratio 

Cxx(T/2)-C X x,ct(T/2) 



7 j w A A V / } AA,Ct V / / ( r\n \ 

flxx = — c^m — • (96) 

i.e. the relative shift of the correlators at T/2. The main unknown here is the coefficient c^- 
Even though it plays a decisive role in the phase diagram of the theory [32], it is difficult to 
obtain in numerical simulations. So far only the ETM collaboration has obtained an estimate 
from their simulations with a twisted mass term [261 El] ■ The data for the pion mass splitting 
together with the LO ChPT prediction gives the rough estimate — 2c2C? ~ (185MeV) 2 at 
a 0.086fm, which translates into | C2 1 ~ (550MeV) 4 . The error, however, is fairly large 
because of the large statistical error in the determination of the neutral pion mass. In 
addition, this value for c 2 was obtained with the tree-level Symanzik improved gauge action 
and the standard Wilson fermion action, and any change in this setup can and probably 
will lead to a different value for c 2 . 9 Nevertheless, for lack of a better estimate we use 
|c 2 | = (500MeV) 4 in the following. 

For the other parameters we use F = 90MeV, a = 0.08fm and a hypercubic lattice with 
Nt = Ni = 24, which corresponds to a box size L = 1.92fm. This implies p ~ 0.75. Even 
though this is slightly large we may still count this as 0(e 2 ) as it should in order to be in 
the GSM* regime. 

Figure [T] shows Rpp and Raa for fx values in the e-regime. For fx = 1.0 we find Rpp = 2.8% 
and it decreases to less than 1% for fx larger than 2. The correction is maximal (less than 
5%) at vanishing jl. However, for fx ~ p we enter the Aoki regime and our formulae are no 
longer valid. The values for Raa are very similar. For instance, Raa — 2.1% at fx = 1.0. 

Figure [l] also shows A a 2(fx), which is, up to the factor p, the numerator in (96). It looks 



very similar to the ratios itself, since the denominator in (96) is of 0(1) and varies only 
mildly for the ft values considered her. So the correction to the correlators is essentially 
A a 2(/i), which happens to be of the order of 10~ 2 . 

The main conclusion we can draw from this exercise is that for our choice of parameters 
the 0(a 2 ) corrections to the correlators are at the few percent level, a comfortably small 
value. 

Using a bigger box improves the epsilon expansion since the expansion parameter 1/(FL) 2 
is smaller. However, a bigger box also leads to larger p values and one easily enters the Aoki 
regime at moderately large volumes. For instance, with Nl = 32 and the other parameters 
unchanged we get L = 2.4 fm and p « 2.4, a value that is certainly not 0(e 2 ). 10 



9 An analysis [5T] of quenched twisted mass lattice data led to a value C2 ~ (300MeV) 4 . 
10 It may seem counterintuitive at first that a change in the volume may bring us into a different regime. 



20 




FIG. 1: Rpp (dotted) and Raa (solid) as a function of fi. Both ratios are smaller than 3.5% for 
jj, > 0.75. The dashed curve represents A a 2 (jl). 

B. Reanalysis of recent lattice data 

In this section we investigate the impact of C2 on the extraction of the continuum low 
energy constants F and E from lattice data. The data is taken from Refs. [221 E2]- It is 
generated with Nf = 2 flavors of dynamical improved NHYP Wilson fermions [53] at a fairly 
small quark mass. From there, a reweighting procedure allows to access even smaller sea 
quark masses [23]- This procedure is exact and does not introduce systematic uncertainties 
but allows to compute correlators at very small quark masses at moderate cost. The lattice 
spacing is a pa 0.115fm from the measurement of the Sommer parameter tq taken to be 
0.49fm[51]. We have two volumes available, one at L — 16a ~ 1.84fm and a larger one with 
L = 24a pa 2.8fm. The former serves mainly as a cross check whereas the latter has sufficient 
size for our NLO formulae to be applicable. Some parameters of the simulation are given in 
TableE 

The theoretical formulae for the pseudo scalar and axial vector correlator both have the 
form constant plus parabola. The coefficient C2 only contributes to the constant term in both 
cases. The curvature itself is rather small at the parameter values simulated, in particular 
compared to the statistical uncertainties, see Fig. [2j Therefore, at a fixed mass, each of 
the two correlators effectively is a constant from which it is difficult to constrain three 
parameters. As already discussed, the theory predicts a particular and relatively strong 



However, increasing the volume requires that we need to decrease the mass in order to stay in the e-regime 
with fixed jl. Hence we have to decrease a as well in order to preserve the relative size between the mass 
and the lattice spacing terms. 
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L/a k ow-pcac A* 



24 



0.128150 0.0019(4) 2.1 

0.128125 0.0024(3) 2.7 

0.128100 0.0030(3) 3.4 

0.128050 0.0044(3) 5.0 
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0.128100 0.0028(11) 0.7 

0.128050 0.0047(9) 1.1 

0.128000 0.0058(7) 1.4 

0.127900 0.0088(5) 2.1 

0.127800 0.0117(3) 2.9 



TABLE I: Parameters of the simulation. L/a is the extend of the box, k the hopping parameters, 
the PC AC quark mass and an approximate values of [i = mT,V, where we use the central value of 



jl dependence of the term multiplied by C2- This gives a handle on the extraction of this 
coefficient. Therefore we simultaneously fit the axial vector and pseudo scalar correlators for 
all available quark masses. From a fit to t e [6, 18] we get £ 1/3 = 249 (4) Me V, F = 88 (3) Me V 
and C2 = 0.02(8)GeV 4 . The data, along with the theoretical curves can be found in Fig. [2} 
Here we used Z|F(2GeV) = 0.90(2) and Z A = 0.99(2) from Ref. j22J. The errors from the 
renormalization factors are not included in the uncertainties of the LECs. The value of C2 is 
compatible with zero within errors and the one sigma band lies within the range of reasonable 
values for a low energy constant. Since the data points are highly correlated, we cannot give 
a good estimate for the quality of the fit; we find x 2 /dof = 0.3(1) without the correlations 
taken into account. We also remark that the results are independent of the fit range once 
tmm/ a > 4. Another concern are the relatively large values of jl. Therefore we repeated the 
analysis leaving the jl ~ 5 data out. We get from the same fit range E 1 / 3 = 250 (4) Me V, 
F = 87(3)MeV and C2 = — 0.01(8)GeV 4 . The differences to the previous values are well 
within the statistical uncertainties. This is encouraging. Even with the additional constant 
the errors of the continuum LECs are reasonably small. 

Is the value we find for C2 large or not? To gauge the impact of this term, we repeat the 
fit by setting C2 = 0. The results are virtually unchanged within errors: X 1 / 3 = 249(4)MeV, 
F = 88(3)MeV. This is very good news. The cut-off effects are so small that they do 
not impact the extraction of the low energy constants beyond the level of the statistical 
uncertainties. 

As a cross check we repeated this analysis on the smaller volume, at the same lattice spac- 
ing and L/a = 16. We obtain X 1 / 3 = 257(4)MeV, F = 83(2)MeV and c 2 = 0.06(14)GeV 4 . 
However, the (uncorrelated) x 2 /dof = 1.3 might indicate that the NLO formulae are no 
longer applicable. These results agree with the findings of Ref. [22J. 
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FIG. 2: Fit of the WChPT predictions to lattice data. All data points within the fit range of 
t/a € [6,18] for the four sea quark masses are included in the combined fit. The axial vector 
correlator is multiplied by a factor 50 for better visibility. 

VI. CONCLUSIONS 

We have shown that the corrections due to the explicit chiral symmetry breaking of 
Wilson fermions are highly suppressed. For typical quark masses these corrections enter at 
either NNLO (GSM regime) or at NLO (GSM* regime). The reason for this suppression can 
be traced back to the fact that the lattice spacing corrections in the chiral effective action 
and the effective operators are either quadratic in a or they come with an additional power 
of either m or p 2 . There is no explicit term with a single power of a only, since such a term 
solely contributes to the additive mass renormalization which is absorbed in the quark mass. 
Hence, the lattice spacing corrections are suppressed in the chiral expansion, similar to the 
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terms in the Gasser-Leutwyler Lagrangian £4. 

In the Aoki regime the modifications are more substantial, affecting the correlators al- 
ready at LO. The main complication in this regime are the zero mode integrals, which are 
no longer the known Bessel functions. 

We tested our formulae against recent lattice data. We found that the additional terms 
which come from the broken chiral symmetry have very little impact on the extracted val- 
ues of F and S whereas the low-energy constant associated with the breaking is hard to 
determine precisely. 

Our results derived here can be generalized in various ways, for example to the case with 
a twisted mass term or to an arbitrary number of flavors. The details of the calculation will 
change, but our various power countings can be carried over with almost no modification. 
Perhaps most interesting from a practical point of view is an extension along the lines of 
Ref. [55], where one considers a mixed setup with some quarks in the e-regime and some 
others in the p-regime. 

However, the main conclusion one can draw is that the effects due to explicit chiral 
symmetry breaking of Wilson fermions in the e-regime are less severe than anticipated before. 
In view of the results of Ref. [22] and the ones presented here, simulations with Wilson 
fermions seem to be a viable alternative to the daunting task of dynamical simulations with 
chiral fermions. 

Note added 

After this paper was completed we received a paper by A. Shindler which also deals with 
Wilson fermions in the e-regime and comes to similar conclusions |56j . 
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APPENDIX A: SOME RESULTS FOR THE EPSILON REGIME 

In this Appendix we summarize formulae which are relevant for the computation of 
correlation functions in the e-regime of chiral perturbation theory. 
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Starting from the leading order continuum chiral Lagrangian of Eq. ([T| and by introducing 



the parametrization of Eq. (25), we can read off the finite- volume scalar propagator for the 
nonzero modes: 



ip;i 



with 



p = 2tt 



no n 



(Al) 



The propagator G(x) satisfies the following properties: 

d 4 xG(x) = 0, 
^6(0) = 0, 
DG(x) = -8(x) ^ 



V 



1 

V' 



(A2) 
(A3) 
(A4) 



UV divergencies, if present, are treated in dimensional regularization. 
We define 03 [57] 

G(0) = A 
T±G(0) s 



T 2 h 



V 

oo 



V 



(A5) 
(A6) 



where ft and &oo are finite dimensionless shape coefficients which depend on the geometry of 
the box. They can be evaluated numerically: for instance, for a symmetric box with L = T 
one has ft = 0.140461 and k 00 = ft/2 (see also [46]). 

In order to obtain time correlators one has to perform integrals over the spatial compo- 
nents of given functions of the propagators G(x). In particular we define [451 EZ] 11 



K^ix - y) = G(x - y)d X/l d Vv G(x - y) - d Xfi G{x - y)d Vv G(x - y) 
+ d XlM d yv H(x-y), 

H{x-y) = / d 4 zG(x- z)G(z-y). 



VJ V 

The integrals that we need for this work are (xq = t): 

t \ T 



(A7) 
(A8) 



d xG(x — y)\ y= o = Thi 



T 



1 

12 



d 3 xK 00 (x - y)\ y =o 



T 

V 



"00, 



dx d yQ J d 3 xH(x-y)\ 



y=0 



T 
V 



hi 



(A9) 

(A10) 
(All) 



11 In the original definition of |45L 157] . K^ v {x — y) contains also contact terms, which we do not consider in 
our computation since we are interested in the correlators for x =/= y. 
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Finally, we recall the SU(Nf) completeness relations which are used for the computation 
of correlation functions. Given the SU(Nf) generators T a , with a = 1, .., iV? — 1 and the 
convention 



one obtains 



Tr [T a T b ] = ^5 ab , 



Tr(T a AT a B) = ]—Tr (AB) + -Tr(A)Tr(B), (A12) 

2N f 2 

Tr(T a A)Tr(T a B) = —Tt(A)Tt(B) + -Tt(AB). (A13) 

2N f 2 



APPENDIX B: SU(2) INTEGRALS 



In the case Nf = 2, the partition function related to the zero-mode integrals in Eq. (54) 
is given by 

Z = [ [dU ]e^ u ° +u b = MM f (Bl) 

JSU{2) f 1 

where I n is the modified Bessel function of the first kind. The normalization 

[dU ] = 1 (B2) 

51/(2) 

has been adopted. Expectation values of arbitrary integer powers of Tr(f/ ) can be obtained 
by computing derivatives of Zq. In particular, for this work we need 

<™o> = ~ = 2|M (B3) 



1 


dZ 




dfx 


1 


d 2 z 


Z 




1 


d 3 Z 


Z 


d[i 3 


1 


d 4 Z Q 


z 





wT/n 3 \ 1 d ' iZo 12 8(3 + /x 2 )/ 2 (2/x) 

; , T . /M 4, 19% 60 24(5 + 2/i 2 )/ 2 (2^) 
((Trt/o) = ^—xt = 16 + -J 3r ,„ v • B6 

Other integrals needed in this work can be related to the previous ones, for instance: 
(TW 2 (TW„) 2 } = ~{TtU ) + 2 (l + \) ((TW ) 2 ) - -<(TW ) 3 } 



/i 

4(15 + 2/i 2 ) 12(10 + 3/i 2 )/ 2 (2/i) 
/i 2 /^i(2//) ' 



(B8) 
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APPENDIX C: OTHER CORRELATORS 



Here we summarize the GSM* results for some other correlators. We start with the 
correlation function of the time component of two vector currents, 



(V a (x)V b (y)} = 6 ab C vv (x-y) 



(CI) 



which we again split into a continuum part and a correction proportional to the lattice 
spacing, 

C vv (x -y) = C VVtCt (x -y) + C vv ^{x - y) . (C2) 
In our notation the leading order vector current in the chiral effective theory reads 



^ct = -i—Tr(T a (U^d,U + Ud^)) 



(C3) 



The continuum contribution at 0(e 6 ) for generic Nf has been calculated before by Hansen 
(see also [19]). After integrating over the spatial coordinates one gets 



ab 



1 T 

-—a v + —koo(3 v 



with 



ay 



(2-U)), 



N 



(C4) 

(C5) 
(C6) 



The function j7o has been defined in Eq. 76 In particular, for Nf = 2 the result explicitly 
reads 



Cy b y(t) 



effy 



»(2/i) \ 



2T 

+ -VTk 00 I 1 - 



1 UM 

fili(2fi)j ' 



(C7) 



The 0(a 2 ) correction in terms of the PCAC mass is given by 



Cyy,a?(t) = -— pA a 2 



(G8) 



where A a 2 is defined in Eq. (95). Comparing this with the result for the AA correlator in 
(93) we observe that the lattice spacing corrections in these two correlators are, up to a sign, 
identical. 

With both the AA and the VV correlator at hand we can trivially obtain the correlation 
functions of right- and left-handed currents. For example, with = [V£ — A^]/2 we find 



C LL (t) = l(Cyy(t) + C AA (t^ 



(C9) 
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and the 0(a 2 ) corrections cancel in the sum on the right hand side, i.e. 
while the continuum part is given by [4"5 

ClL,ctv) 



(CIO) 



T 
2T 



1 & T i. 



V h(2fj) 1 \T 



(Cll) 



The same result can be obtained by a direct calculation of the correlator, of course. 
Finally, the scalar correlator 

(S a (x)S b (y)) = 5 ab C ss (x - y) (C12) 

vanishes identically in the chiral effective theory for Nf = 2, as one can check either by 
explicit calculation or by using G-parity. 
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